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ABSTRACT 


It  is  demonstrated  mathematically  that  any  desired  form  function 
can  be  approached  by  the  burning  of  increments  or  groups  of  dual 
composition  propellant  grains  of  neutral  geometry  (constant  surface 
during  burning).  This  scheduled  burning  is  a.chieved  by  varying  the 
thickness  of  each  of  the  outer  layers  of  propellant,  with  the  outer 
layers  having  the  lower  burning  rate. 

It  is  shown  that  the  pertinent  parameters  are  the  total  burning 
distance  of  the  propellant  system,  the  burning  distances  of  grains  in 
each  increment  of  propellant,  the  nximber  of  groups  or  increments, 
the  mass  of  each  group,  and  the  number  of  groups  which  are  burning 
simultaneously.  The  relationship  of  these  parameters  to  the  desired 
form  function  is  expressed  as  a  series  of  equations,  each  of  which  is 
an  incremental  form  function.  The  sum  of  the  series  approaches,  to 
any  desired  degree  of  accuracy,  the  specified  form  function. 
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ABSTRACT 


It  is  demonstrated  mathematically  that  any  desired  form  function 
can  be  approached  by  the  burning  of  increments  or  groups  of  dixal 
composition  propellant  grains  of  neutral  geometry  (constant  surface 
during  burning).  This  schedioled  burning  is  achieved  by  varying  the 
thickness  of  each  of  the  outer  layers  of  propellant,  with  the  outer 
layers  having  the  lower  burning  rate. 

It  is  shown  that  the  pertinent  parameters  are  the  total  burning 
distance  of  the  propellant  system,  the  burning  distances  of  grains  in 
each  increment  of  propellant,  the  number  of  groups  or  increments, 
the  mass  of  each  group,  and  the  number  of  groups  which  are  burning 
simultaneously.  The  relationship  of  these  parameters  to  the  desired 
form  function  is  expressed  as  a  series  of  equations,  each  of  which  is 
an  incremental  form  function.  The  sum  of  the  series  approaches,  to 
any  desired  degree  of  accuracy,  the  specified  form  function. 
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INTRODUCTION 


Propellant  devices  generally  produce  a  driving  gas  through 
chennical  reaction.  The  number  of  moles  of  gas  produced  depends 
on  propellant  composition,  while  the  mass  consumption  rate  is, 
additionally,  a  function  of  propellant  grain  geometry  and  pressure. 

At  constant  pressure,  constant  linear  burning  rate  makes  the  mass 
consumption  rate  a  function  of  total  burning -surface  area.  The  re¬ 
lation  between  mass-fraction  burned  and  distance -fraction  burned 
is  unique  for  every  grain  geometry.  This  relation  is  convention¬ 
ally  and  appropriately  called  the  "form”  fvuiction.  Thus,  each 
grain  configuration  will  exhibit  a  specific  gas  generation  schedule. 

It  has  been  demonstrated*  that  most  form  functions  can  be 
satisfactorily  represented  by  a  qua-dratic  equation  of  the  form 

N/C  :=  ko  +  kiZ  +  kgZ^  (1) 

where  N  =  mass  of  propellant  cons  ^med 
C  =  total  mass  of  propellant 
Z  =  distance  fraction  consumed 
k  =  constant. 

As  used  here,  the  burning  distance  is  defined  as  the  shortest 
distance  normal  to  the  burning  surface  that  a  grain  recedes  until 
it  either  loses  its  structural  integrity  or  is  completely  consumed  - 
whichever  comes  first.  The  form  function  for  a  neutral  geometry 
grain  (such  as  a  sheet  propellant)  is 

N/C  -  Z  (2) 

It  will  be  demonstrated  that  any  generalized  form  function,  g(Z), 
based  on  a  burning  distance  W,  can  be  approached  by  the  incremental 
addition  of  several  form  fxinctions  of  the  form  N/G  =  Z,  based  on 


*C,  F.  Curtiss  and  J.  F.  Wrench,  "Interior  Ballistics,"  Carnegie 
Institute  of  Washington,  15  July  19^5 
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groups  of  neutral  geometry,  sheet-type  grains  of  burning  distance  w, 
where  w  <  W.  Neutral  geometry  grains  are  defined  as  having  a  con¬ 
stant  burning  surface  during  consumption. 


THE  DISTRIBUTED  PROPULSION  SYSTEM 

The  total  propellant  charge  is  composed  of  n  groups  (n  is  an 
integer)  of  neutral  geometry  grains.  The  mass  per  group  is  Cj, 
where  i  runs  from  i  =  1  to  i  =  n.  All  grains  in  any  one  group  will 
commence  burning  simultaneously  when  the  burning  distance  of  the 
grains  in  the  preceding  group  has  receded  to  sorr^e  fraction  w/a. 

The  parameter,  a,  is  defined  as  the  reciprocal  of  the  distance  frac¬ 
tion  of  the  first  group  burned  at  the  time  group  two  begins  to  burn. 
Furthermore,  the  last  (a)  groups  will  be  completely  consumed  si¬ 
multaneously.  Thus,  the  over-all  distance  burned,  W  (the  total 
distance  burned  during  consumption  of  the  entire  charge),  is  nw/a. 

Figure  1  is  a  graphical  illustration  of  the  mathematical  model 
for  the  case  of  nine  groups  of  grains  (n  =  9)  with  each  group  com¬ 
mencing  burning  when  1/3  of  the  preceding  group's  burning  distance 
has  been  consumed  (a  =  3).  Here,  the  abscissa  represents  the  burn¬ 
ing  distance  of  grains  in  each  group.  When  a  group  commences 
burning  (for  example,  group  6),  some  groups  may  have  been  com¬ 
pletely  consumed  (groups  1,2,  and  3),  other  groups  partly  con¬ 
sumed  (group  4  is  2/3  consumed  and  group  5  is  1/3  consumed),  and 
other  groups  will  not  have  begun  burning  (groups  7,  8,  and  9). 

Note  that  the  over-all  burning  distance,  W,  is  equal  to  n/a  (or 
three)  times  the  burning  distance,  w,  of  the  first  group.  Note  also 
in  Figure  1  that,  for  the  last  (a)  groups  to  be  consumed  simviltaneously, 
the  burning  distances  of  the  latter  (a  -  1)  groups  (groups  7  and  9)  must 
be  less  than  the  burning  distance,  w,  of  the  first  group.  The  reason 
for  using  neutral  geometry  grains  is  that  their  form  function  is  linear 
and,  thus,  when  individual  groups  of  grains  are  burned  as  described 
above,  the  over-all  form  function  can  be  expressed  as  a  series  of 
linear  segments. 
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Graphical  Illustration  of  the  Programmed  Burning  Model 


Mathematical  Description  of  the  Model 


Consider  a  propellant  system  with  n  groups  of  scheduled  grains, 
such  that  the  ith  group  is  commencing  burning.  If  i  >  a,  the  first 
(i  -  a)  groups  will  have  been  completely  consumed.  Any  other  groups 
preceding  group  i  will  have  been  partly  consumed.  All  other  groups, 
including  group  i,  will  not  have  commenced  burning.  Figure  1  illus¬ 
trates  the  case  where  a  =  3,  and  i  is  some  intermediate  group  (e.g.  , 
i  =  6). 

The  total  mass  fraction  of  core  propellant  consumed  as  group  i 
commences  burning  can  be  expressed  in  summation  form  as 


n+1  -i 

L  C.  +  S  [(a 

^  k=l 


k)/a-^  ^i-a+k 


(3) 


i+a-(n+2) 

+  L  [(i  -  n  -  q  +  a  -  l)/(a  -  q)J  C 

q=l 


n+q-a+1 


J 


where  4  =  the  form  function  for  the  distributed  propulsion  system. 


The  subscripts  on  the  C's  refer  to  the  integral  number  of  the 
group.  The  first  siimmation  term  represents  the  first  (i  -  a)  groups 
that  are  entirely  consumed.  Restrictions  on  i  and  j  are  that  i  must 
be  greater  than  (a),  and  j  greater  than  one.  The  second  summation 
term  represents  those  (n+1  -  i)  groups  preceding  group  i  that  are 
partially  consumed  and  whose  burning  distance  fractions  consumed 
can  be  expressed  in  integral  multiples  of  1/a,  The  restriction  on 
this  second  term  is  that  k  must  be  greater  than  (a  -  i).  The  last 
summation  term  represents  those  groups  preceding  group  i  that  are 
partially  consumed  and  whose  burning -distance  fractions  consumed 
cannot  be  generally  expressed  in  terms  of  integral  multiples  of  1/a. 
Three  restrictions  on  this  term  are  that  (a)  be  greater  than  one,  and 
that  q  be  greater  than  (a  -  n  -  1),  and  less  than  (a). 


Statement  of  the  Problem 


A  solid  propellant  device  will  have  associated  with  it  a  specified 
propellant  form  function,  g(Z),  and  a  burning  distance  W.  The  prob¬ 
lem,  in  terms  of  a  distributed  propulsion  system,  can  be  stated  as 
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follows:  based  on  a  burning  distance  of  W,  using  n  groups  of  neutral 
geometry,  sheet  type  grains  such  that  a  maximum  of  (a)  groups  burn 
at  a  time,  generate  form  function,  g(Z);  then,  determine  the  mass 
and  burning  distance  for  each  group  based  on  a  core  propellant  com¬ 
position  which  is  identical  with  that  required  in  the  propellant  system. 


Solution  of  the  Problem 

Since  both  n  and  (a)  have  been  specified,  the  burning  distance 
w,  for  grains  in  all  groups  except  the  latter  (a  -  1)  groups  can  be 
determined  from  the  relation 

w  =  (a/n)W  (4) 

This  relation  becomes  more  evident  if,  in  Figure  1,  a  comparison 
is  made  between  the  over-all  burning  distance  W  and  the  individual 
burning  distance  w.  Also,  from  Figure  1,  the  burning  distance  for 
grains  in  the  latter  (a  -  1)  groups  are,  successively, 

(a  -  1)W  ,  (a  -  2)W . [a  -  (a  -  l)3w  (5) 

n  n  n 

The  total  mass  fraction  burnt  can  be  calculated  from  Equation  3  at 
n  values  of  i,  namely  1,  2,  3,  .  .  .  .  , i , .  .  .  ,  n.  These  values  corre¬ 
spond  to  n  definite  values  of  Z,  by  the  equation  Z  -  (i  -  l)/n,  namely 

0/n,  1/n,  2/n . .  {i  -  l)/n,  (n  -  l)/n  (6) 

From  Equation  3,  the  n  values  of  mass  fraction  consumed  are 


^1,  ^2,  . .  ^i, . .  ^n.  (7) 

For  each  value  of  Z,  a  corresponding  value  of  g(Z),  the  specified 
mass  fraction  consumed,  can  be  generated: 

g{0/n),  g(l/n) ,  g(2/n) ,  .  .  .  ,  g  [{i  -  l)/n]  ,  .  . .  ,  g  Qn  -  l)/n]  (8) 


5 


Upon  setting  the  specified  and  the  scheduled  mass  fractions 
equal  to  each  other  for  correspondence  values  of  Z  and  i,  a  set  of 
n  independent  equations  with  n  unknowns  is  generated. 


(n-l)/n 

n-1 

g(Z)  = 

Z=0 

i=0 

g(i/n)  = 


4+1 


(9) 


which  can  also  be  expressed  as 


n-l| 


i=0 


g{i/n)=  Cj  +  0!i+i^2  ^2  “i+l.i+1  ^^i+l 


(10) 


+  •  •  •  •  +  Oi.  ,  I  ^n 
1+1, n  ^ 


where  the  OL  terms  are  the  coefficients  of  the  C's  in  Equation  3, 
An  additional  equation  applies  after  complete  consumption  of  all 
groups  of  grains; 


g(n/n  =  1)  =  i.  I  ^ 
C  l^pHl 


(11) 


where  p  denotes  the  group  number.  This  additional  equation  is 
necessary  for  solution  of  the  set,  since  the  solution  to  the  first 
equation  in  the  set,  g{0/n)  =  is  trivial. 


The  general  solution  of  the  set  is 


a.  O' 

i+1,1  i+1,2 


O'. 


i+l,n 


(12) 
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where  D  is  the  determinant  of  the  coefficients. 

Graphical  Description  of  the  Model 

Figure  2  illustrates  how  the  incremental  addition  of  the  form 
functions  of  five  groups  of  grains  results  in  some  selected  form 
function.  This  figure  represents  a  comparison  of  the  over-all  form 
fvinction,  where  a  =  2,  with  individual  form  functions  for  each  group. 
The  mass  fraction  consumed  at  any  point  on  the  over-all  form  func¬ 
tion  equals  the  summation  of  mass  fractions  of  all  the  groups  either 
totally  or  partially  consumed  at  that  point.  Note  that  the  total  mass 
fraction  consumed  at  any  jointure  point  (points  marked  X)  depends 
on  the  mass  fraction  consumed  at  the  preceding  jointure  point,  plus 
the  contributions  of  the  (a)  groups  (here,  a  =  2)  burning  simultane¬ 
ously  in  the  preceding  interval  of  distance  fraction  consiimed.  This 
interval  is  w/a  (or,  in  Figure  2,  w/2)  in  length. 

Excluding  the  latter  (a 1)  groups,  it  is  generally  true  that 
the  slope  of  the  form  function  for  each  individual  group  is  propor¬ 
tional  to  the  mass  of  that  group.  Thus,  groups  three  and  four, 
which  have  identical  slopes  in  Figure  2,  can  be  expected  to  have  the 
largest  masses  of  the  first  four  groups. 

Illustration  of  the  Use  of  the  Distributed  Propulsion  System 

2 

Afssume  that  it  is  desired  to  duplicate  the  form  function,  g(Z)  =  Z  , 
for  a  large  gun  system  consisting  of.  4,  0  pounds  of  propellant  divided 
into  1334  grains  of  0,  2-inch  buraiag  distance,  with  a  distributed  pro¬ 
pulsion  systan  consisting  of  10  groups  of  sheet  type  grains  of  identical 
composition.  Let  the  maximum  number  of  groups  burning  at  one  time 
vary  from  one  to  five.  The  problem  is  to  calculate  the  mass  and  burn¬ 
ing  distance  of  each  group  for  each  value  of  (a)  (a  varies  from  one  to 
five,  and  n  equals  10).  Solution:  for  n  =  10,  a  =  1,  Equation  3  becomes 

C  $  =  C  4  'l  4  '  e'  Ci04q 

j  =  l  q=l  ^  ^  ' 
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.  Incremental  Addition  of  Form  Functions 


Corresponding  values  for  Z  and  g(Z)  when  group  i  commences 
burning  are  presented  in  Table  I.  Note  that  Z  is  determined  from 
the  relation,  Z  =  {i-l)n,  and  is  independent  of  the  value  of  (a). 


TABLE  I,  Values  of  the  Form  Function,  g(Z)  =  Z^ 


Group  i 

c 

1 

B 

B 

<=4 

c 

5 

B 

c 

i 

G 

8 

c 

9 

B 

Z 

0 

B 

B 

B 

B 

.6 

B 

.8 

B 

g(Z)  =  Z^ 

0 

•  01 

B 

.09 

B 

B 

.  36 

.49 

.64 

.81 

An  additional  equation  is  applicable  after  complete  consumption: 


N 

C 


1.0 


Cl  + 


+  G 


10* 


Upon  equating  the  g(Z)'s  to  the  <i?'s  for  cor  responding  values  of 
Z  and  i,  and  solving  the  set  of  simultaneous  equations  so  formed, 
the  following  group  masses  (Table  II)  E-re  obtained. 


TABLE  II,  Group  Masses 


Group 

c 

2 

c 

3 

c 

4 

c 

5 

c. 

6 

C 

•n 

i 

c 

8 

c 

9 

c 

10 

Mass  (lb) 

.04 

.  12 

.  20 

.28 

.36 

.44 

.52 

.  60 

.  68 

.76 

Burning  distances  can  be  calculated  from  Equations  4  and  5, 
Table  III  offers  a  comparison  of  the  masses  and  burning  distances  of 
each  group  for  various  values  of  (a).  Note  that  the  latter  (a  -  1) 
groups  have  smaller  burning  distances  than  the  other  groups. 
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TABLE  III,  Mass  (lb)  and  Burning  Distance  (in.)  of  Each  Group 
for  a  =  1  to  a  =  5 


Group 

Value 

■=1 

^2 

S 

^4 

S 

^6 

^8 

S 

=  10 

M 

.04 

.  20 

.  28 

.  36 

.44 

.52 

.  60 

.  68 

.76 

1  W 

.02 

.02 

.02 

.02 

.02 

.02 

.  02 

.02 

.  02 

.02 

,08 

.  16 

.24 

.32 

.40 

.48 

.  56 

.64 

.72 

.40 

2 

.04 

.0^ 

.04 

.04 

.04 

.04 

.04 

.04 

.04 

.02 

.  12 

.24 

.24 

.36 

.48 

.48 

,  60 

.72 

.48 

.  28 

3 

.06 

.  06 

.06 

.06 

.06 

.06 

.06 

.06 

.04 

.02 

.  16 

.32 

.32 

.32 

.48 

.64 

.64 

.48 

.40 

.24 

4 

.08 

.08 

.08 

.08 

.08 

.08 

.  08 

.  06 

.04 

.02 

.20 

.40 

.40 

.40 

.40 

.60 

.  64 

.  48 

.  32 

.  16 

5 

.  10 

.  10 

.  10 

.  10 

.  10 

.  10 

CO 

o 

• 

.06 

.04 

.02 

M  -  Mass  (lb) 
W  -  Web  (in) 


Limitations  on  a  and  n 

Certain  limitations  exist  on  the  values  of  n  and  (a).  From 
Equation  4,  it  is  evident  that  w  decreases  as  n  increases.  Thus,  the 
limiting  (smallest  practical)  value  of  w  determines  the  largest  value 
of  n.  The  smallest  value  of  n  is  determined  from  the  desired  degree 
of  approach  of  the  scheduled  form  function  to  the  specified  form  func¬ 
tion,  As  n  increases,  the  distributed  form  function  approaches  the 
specified  form  function  more  closely. 

Other  factors  that  may  impose  limitations  on  n  and  (a)  are  grain 
formulating  and  forming  variables  encountered  in  manufacturing  and 
certain  aspects  of  the  propellant  system,  such  as  chamber  geometry, 
ignition,  and  loading  density.  These  factors  will  vary  with  the  propel¬ 
lant  used  and  the  individual  test  setup. 
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Grain  Prepa.ration  Techniques 

It  is  envisioned  that  an  outer  layer  of  propellant  with  a  much 
slower  burning  rate  than  that  of  the  core  propellant,  completely 
covers  all  available  grain  surface  area.  The  thickness  of  the 
outer  propellant  for  a  specific  grain  group  will  depend  on  the 
scheduling  time  to  the  commencement  of  burning  for  that  group. 
Hereafter,  this  outer  propellant  layer  will  be  referred  to  as  the 
scheduling  propellant. 

Two  grain-forming  techniques  are  being  employed  in  current 
experimental  work.  The  first  technique  utilizes  preformed  sheets 
of  two  double  base  propellants  having  widely  different  burning  rates, 

A  "sandwich”  grain  is  formed,  with  the  faster  burning  sheet  lying 
between  two  sheets  of  the  slower  burning  scheduling  propellant. 
Uncured  core  propellant  is  employed  as  the  adhesive.  Thus,  succes¬ 
sive  groups  of  grains  will  be  characterized  by  increasing  thicknesses 
of  the  scheduling  propellant.  As  the  thickness  of  the  outer  layer  in¬ 
creases,  the  burning  time  to  reach  the  core  propellant  increases  at 
constant  pressure. 

The  second  technique  involves  the  application  of  a  nitrocellulose 
lacquer  to  a  sheet  of  core  propellant.  The  basic  variables  to  be 
considered  here  are  the  relative  percentages  of  added  resins,  plas¬ 
ticizers,  active  and/or  inactive  solvents,  and  diluents  in  the  lacquer, 
application  techniques,  and  nitocellulose  layer  thickness. 

It  has  been  estimated  that  when  the  burning  rates  at  constant 
pressure  of  the  two  propellants  differ  by  an  order  of  10  :  1  (with  the 
outer  layer  obviously  having  the  slower  burning  rate),  the  outer,  or 
scheduled,  propellant  composition  is  only  seven  percent  by  weight  of 
the  propellant  system.  Thus,  properties  of  the  gas  formed  are  largely 
determined  from  the  properties  of  the  core  propellant  alone. 

One  criterion  for  the  dual  compositions  in  any  system  operating 
over  a  range  of  pressures  is  that  the  ratio  of  their  burning  rates  be 
constant  and  independent  of  pressure.  Interpreted  mathematically, 
this  criterion  states  that  the  pressure  exponents  in  the  burning  rate 
expression  for  each  of  the  two  propellant  compositions  should  be 
equal. 

The  masses  C  ,  etc,  ,  represent  masses  of  the  core  propel¬ 

lant  only.  The  method  of  calculation  of  the  mass  of  scheduling 


11 


propellant  in  each  group,  for  a  specified  form  function,  will  depend 
on  the  properties  of  the  scheduling  propellant  -  specifically,  on  the 
burning  rate-pressure  relationship  and  the  density  of  the  propellant. 

Results  on  the  burning  behavior  of  the  dual  composition  grains  in 
a  closed  bomb  and  a  vented  vessel  apparatus  will  indicate  the  effective¬ 
ness  of  these  techniques  in  providing  scheduled  burning  in  accordance 
with  the  model.  These  results  will  basically  consist  of  pressure-time 
data  and  dP/dt  vs  time  data. 


CONCLUSIONS 

1.  Any  desired  form  fiinction  can  be  approached  by  the  scheduled 
burning  of  a  number,  n,  of  increments  of  dual  composition,  neutral 
geometry,  solid  propellant  grains.  The  form  function  so  achieved 

is  dependent  on  the  mass  of  the  inner  (core)  propellant  for  each  in¬ 
crement,  the  thickness  for  each  increment  of  the  outer  (scheduling) 
propellant  covering  grains  of  uniform  size  in  that  increment,  and  the 
ratio  of  burning  rates  of  the  two  propellants. 

2,  The  degree  of  approach  of  the  scheduled  form  function  to  the 
desired  form  function  increases  as  n,  the  number  of  increments 
(groups)  is  increased.  As  the  number  of  increments  burning  simixl- 
taneously  (a)  is  increased,  the  burning  distances  of  the  core  propel¬ 
lant  becomes  more  nonuniform  from  increment  to  increment.  In 
general,  the  latter  (a  -  1)  increments  have  smaller  burning  distances 
than  the  other  increments. 


RECOMMENDATIONS 

It  is  recommended  that  tests  be  carried  out  to  test  the  basic  model. 
It  is  expected  that  the  results  of  such  tests  will  reveal  which  of  the  two- 
grain  preparation  techniques  -  laminating  propellant  layers  or  applying 
nitrocellxxlose  lacquer  films  -  is  better  from  a  practical  standpoint. 
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